Quantum transitions induced by the third cumulant of current fluctuations 
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We investigate the transitions induced by non-Gaussian external fluctuations on a small quantum 
system. The rates for the transitions between the energy states are calculated using the real-time 
Keldysh formalism for the density matrix evolution. We detail the effects of the third cumulant 
of current fluctuations coupled to a quantum system with a discrete level spectrum and propose a 
setup for detecting the frequency-dependent third cumulant through the transitions it induces. We 
especially discuss a scheme where the fluctuations are coupled to a Josephson flux qubit. 



The study of fluctuations has been in the center of 
interest in physics for decades. The relevance of noise 
and fluctuations is underlined by the fundamental re- 
lation between fluctuations and dissipation in physical 
systems. One very concrete example of fluctuations is 
the current noise in electric circuits. At equilibrium, it 
obeys the fluctuation-dissipation theorem which relates 
the magnitude of fluctuations to the temperature and the 
impedance of the circuit. For a quantum system with a 
finite number of levels interacting with an environment, 
the magnitude of these fluctuations in the environment 
then determines the steady state of the system, along 
with the rate with which this steady state is approached. 

During the past decade, the theory of electric fluctu- 
ations in mesoscopic systems has been significantly de- 
veloped to characterize them also out of equilibrium^^ 
where a finite average current leads to shot noise. The 
study yields information about the microscopic physical 
phenomena inside electric conductors and the effects of 
the electromagnetic environment on mesoscopic circuits. 
In large wires the current statistics is Gaussian and fully 
characterized by the average current and the noise power. 
The experimental development in manufacturing smaller 
circuits has enabled the study of the non-Gaussian char- 
acter of fluctuations in mesoscopic samples In princi- 
ple, the knowledge of these fluctuations allows for an im- 
proved characterization of the conductors^ or the study 
of the effect their non-Gaussian character causes on other 
mesoscopic systems ^LLSi 

With a nonvanishing average current, the probability 
distribution of current fluctuations no longer needs to be 
symmetric around the average current. In particular, the 
third cumulant of fluctuations describing the skewness 
of the current distribution may be finite. It is also the 
lowest cumulant indicating a non-Gaussian distribution. 
Despite the strong theoretical effort describing the nature 
of the higher-order cumulantsji measuring even the third 
cumulant with conventional techniques has turned out to 
be difficult and so far its only measurements exist for the 
case of a tunnel junction^ The attention is thus turning 
towards using other mesoscopic systems as fluctuation 
detectors, 5 i 6 i 9 i 10 i n 

In this Paper we analyze the transitions caused by ex- 
ternal fluctuations on a probe quantum system. First, we 
present a formula correcting the Golden Rule transition 
rates by taking into account the next order effects that 



are dependent on the third cumulant. This is essential in 
developing generic methods for detecting non-Gaussian 
fluctuations. We can establish conditions imposed to 
suitable probes of third-cumulant induced excitations. 
Although we concentrate on current fluctuations, our 
general analysis is independent of the physical system 
as long as the fluctuations are linearly coupled to the 
probe system. To demonstrate the results, we consider 
a quantum two-state system (qubit) as a probe candi- 
date and propose a setup for measuring the effects of the 
frequency-dependent third cumulant of current fluctua- 
tions by a Josephson flux qubit>i2*i» This can be viewed 
as a generalization of using qubits as spectrometers of 
the quantum noise power^i a method which has already 
been experimentally demonstrated^* 
Our starting point is the Hamiltonian 

H = ifext + H s + flint, (1) 

where fl ex t and fl s describe the environment where the 
current fluctuates and the quantum system we use as a 
probe for the fluctuations, respectively, and fl; nt is the in- 
teraction Hamiltonian between the environment and the 
probe. Motivated by the case of a current-biased Joseph- 
son junction and the magnetic interaction between two 
circuits considered below, we study the bilinear coupling 
of the form H lnt = g SIcf>. Here 51 is the current fluctua- 
tion operator acting on the environment, <fi is an operator 
acting on the probe system and g is the coupling constant 
of the interaction. We assume that the quantum system 
is described by a set of energy eigenstates {\n}} and the 
average current effect g (I)4> is included in H s . Treating 
flint as a perturbation, the Fermi Golden Rule predicts 
the transition rate T^\ n , = ^§-\(j) n n>\ 2 S S i( En ~® n ' ) be- 
tween the eigenstates of the probe systemiAi The matrix 
element is defined as <p nn > — (n\4>\n') and the noise power 
S 5i (uj) = i J^e iwt (<5/(t)<5I(0))cZi. The lowest-order 

estimate r^ 2 ) is thus proportional to the second cumu- 
lant of current fluctuations. The correlator in the above 
expression is calculated with respect to the environment 
Hamiltonian H oxt as if the probe system did not exist. 
Below, we correct the transition rate by calculating 
the next order contribution T^ 3 \ depending on the third 
cumulant. 

We solve the density matrix evolution using the real- 
time Keldysh method, as outlined in Refs. lilmri which is 
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a natural formalism for studying a small subsystem in a 
larger environment. We are interested in the dynamics of 
the probe system in particular, so we study the reduced 
density operator p(t) = Tr cxt ptot{t) where p to t is the den- 
sity operator for the system and the environment. The 
trace goes over a complete set of environment states. The 
idea is to solve the temporal evolution of a diagonal ele- 
ment of the reduced density matrix p n > n i (t) = (n'\p(t)\n') 
with the initial condition p nn (to) = 1 (n ^ n'). In 
the long-time limit p n > n '(t) is proportional to the total 
evolution time t — to, the coefficient being the transi- 
tion rate Y n ^ n i. We calculate the rates between well- 
specified states of the reduced system. Therefore, with- 
out loss of generality, we use an initial state of the form 
Ptot(io) = Poxt(to) ®p(<o)> where /9 ex t(*o) describes the 
initial state of the environment. 



FIG. 1: Six diagrams contributing to r^ n , . Each diagram 
represents temporal evolution from the initial state to 
the final state \n')(n'\. The dots represent interaction vertices 
g 51(f) and the horizontal lines forward and backward propa- 
gation in time. 

As the lowest-order contribution to T n ^ n i is the well- 

(2) 

known Golden- Rule result T y n ^ n , , we concentrate on the 

(3) 

next order contribution T n ^ n , . The total rate is then 
given by T n ^ n , = T^l n , +T^l n , . Treating H int = g 81<j> 
as a perturbation and using the graphical rules derived 
in Ref. [Til we find six different diagrams contributing to 
if^,, see Fig. □ 

To present the result in a compact way we define a 
correlator 
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d{t 3 -h) / d{t 2 -h)x 



x gtuite-to-H^te-ti^/f [5I(ti)5Ifo)]5I(ta)), (2) 

where T denotes the anti-time-ordcring operator. The 
time-dependent correlator is calculated with respect to 
the free external Hamiltonian H ext with the density oper- 
ator yOext(^o)- We assume H cxt to be independent of time 
and Pext(io) to describe a stationary state with respect 
to -ffext- Our results can be also stated with the help of 
the Fourier transform of (8I(t3)T[6I(t2)SI(ti)]) , which is 
the complex conjugate of the previous correlator, so it is 
a matter of choice which one to use. Our definition of 
the frequency-dependent third cumulant J2J) differs from 
the one studied in Ref. 0, which consists of the sum of 
all possible Keldysh orderings. Whereas that definition 



is relevant in studying the evolution of the off-diagonal 
density matrix dynamics, transition rates cannot be ob- 
tained from that form. If the state of the environment 
is invariant under time reversal as usually in equilibrium 
at low magnetic fields, the correlator vanishes. 

Evaluating and summing the different contributions 
shown in Fig. ^ we obtain the result 
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The summation is extended over all the eigenstates of H s 
and rj denotes a positive infinitesimal quantity. With the 
help of the identity x _ x \± iri = P^^T^S(x-x ), where 
P stands for a principle value integral, we can write © 
in the form 
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Even without the knowledge of 8 3 I(uii, UJ2), the general 
results llol l[l contain some information about the require- 
ments made for the meter designed to detect the third- 
cumulant effects. The structure of the product of the 
matrix elements 0„ irii restricts the possible physical re- 
alizations used in detecting the transitions induced by 
the third cumulant. Generally the operator <f> should ei- 
ther couple several states of the system, or both matrix 
elements 4> n ,n and <j) n ,n+x should be finite. 

Next we turn to study the case where the probe sys- 
tem is a qubit. The system Hamiltonian can be written 
as H s = —\B z <Jz — \B x o x and the interaction term as 
^int — 9 81 a 2 . The system Hamiltonian has the eigen- 
states \Ex) = a\ T}+/3| |), \E ) = -(3\ T) + a| |) and the 
eigenenergies E x = \^/Bj + B 2 z , E = -\^/Bl+B 2 z . 
The coefficients can be parametrized as a — cos^ and 
(3 = sin§, where <j> = arctan(|p). We denote the en- 
ergy difference between the two eigenstates as AE = 
y/B 2 + B 2 . Using the above conventions and the general 
result 10}, we can express the corrections to the transition 
rates as 
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The function F(uj) contains the information about the 
third cumulant and is defined as 
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Comparing the result (JSJ) with the Golden Rule rates 
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one notices that the function F(oj) plays a similar role in 
r( 3 ) as the noise power m 1 y ' . 

Supposing we can control the effective magnetic fields 
B x and B z , we can optimize the parameters a and (3 
to produce the maximum effect from I^ 3 ). The abso- 
lute value of the expression (a /3) 2 (a 2 — f3 2 ) is maximized 
by choosing a = 0.89 and j3 = 0.46 or vice versa, i.e., 
B x = 1AB Z or B z = \AB X . By changing the magnitude 
of AE = ^B 2 + B 2 Z but keeping B x /B z fixed, one can 
probe F(uS) as a function of frequency. 

A physical qubit always has some intrinsic noise mech- 
anism, in solid-state realizations produced by the elec- 
tromagnetic environment, which cannot be neglected (we 
consider the external fluctuation circuit as an additional 
environment). To be measurable, the external current 
fluctuation effects have to be significant compared to 
transitions due to the intrinsic noise. 
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FIG. 2: Josephson flux qubit inductively coupled to an exter- 
nal circuit producing current fluctuations. The effective mag- 
netic fields B Z ,B X can be tuned by controlling fluxes $i,$2 
to maximize the effects of . 

A possible physical realization for the system consid- 
ered above is a Josephson flux qubi t 12 ' 13 coupled induc- 
tively to the external circuit, see Fig. The interaction 
Hamiltonian is of the form H lnt = Y T * 51 a z , where M 
is the mutual inductance between the qubit and the ex- 
ternal circuit, A<f> is the flux difference between the two 
states of the flux qubit and L q t is the inductance of the 
qubit. We choose (SI) = 0, since the effects of the finite 
average external current can be included in redefining B z 
or eliminated by a flux control. The effective magnetic 
fields can be controlled by external fluxes through the 
loops so the qubit can be biased to the optimal point for 
detecting T^ 3 ). The transition rates follow from Eqs. 
and (JZJ) after the identification g — • We- assume 
that the energy gap to the higher states is large com- 
pared to any other energy scales in the system, allowing 
us to make the two-state approximation and to neglect 
the effective interaction terms nonlinear in 51. 

Let us estimate I^ 3 ) in a flux qubit for a specific setup. 
Suppose that the external circuit consists of a scatterer 



with resistance R and loop inductance L. We assume 
that the third cumulant of current fluctuations in the 
scatterer is frequency independent in the frequency scale 
of the circuit, lol = R/L. This is generally the case 
provided that the voltage eV over and the Thouless en- 
ergy Et of the scatterer, defined as the inverse time of 
flight through it, satisfy eV, E T > huj L . 17 - 18 Then the 
frequency dependence of the correlator arises solely 
from the classical effect of the inductance L modifying 
the noise. In this limit Eq. @ can be approximated by 
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where / is the average current in the circuit and F3 is 
a scatterer-specific proportionality constant ("Fano fac- 
tor") between the third cumulant and the current. 

In deriving JSJ, we assumed L q bl 2 b <C LI 2 where I q i, is 
the current in the qubit, allowing us to neglect the back- 
action of the qubit on these fluctuations. This leads to 
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r' 3 ' given by 
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where A = 32TrF 3 e 2 Ig 3 (a(3) 2 (j3 2 — a 2 ). The noise power 
for the setup can be written as 
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where F2 is the Fano factor for the second cumulant and 
0(x) is the Heaviside step function. This formula includes 
the quantum fluctuations (last term) and is valid for our 
case provided that the temperature T is low, UbT <C AE. 
In the limit eV > AE we get from Eqs. Q), 10 and JTUJ) 
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with g = ( A ^^ e ) • For the optimal parameters a and (3 

mentioned above, a 2 — 1 = 0.58. The phase difference 
of the two flux states can be of order <&o/4 = h/8e, so 

we may estimate g sa ( sz/f ) ' Consequently it can be 
made of order unity or greater by an efficient inductive 
coupling and a large external inductance L. The factor 
F3/F2 depends solely on the nature of noise produced 
by the scatterer^ For realistic parameters lol/{2tt) 
1(1 GHz and AE/h = 1 GHz the last factor is about 2.5%. 
Optimizing the setup one could expect a relative effect 

Y^aj up to roughly 10 %, which shows that the third 

cumulant effect can be significant. 

Now suppose that the intrinsic relaxation of the qubit 
is caused by an independent zero-averaged fluctuating 
Gaussian field. Then the second-order rate should be 
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replaced by the sum of rates caused by the field and 
the external circuit, the third-order rate remaining un- 
changed. In the case of a zero-temperature environment, 
this intrinsic relaxation rate T^t can be quantified by 
the (J-factor, r int = AE/hQ. In this case, its ratio to 
the rate T^ 2 ) is 
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Here Rq = h/e 2 . 

One possibility to detect I^ 3 ) is to let the qubit reach 
the stationary state and then determine the probabilities 
P Eo and P El = 1 - Pe of the states \E ) and \Ex). This 
can be achieved by repeated measurements of the qubit. 
From detailed balance we get 
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where 7Q = 7 >f + 7 £ and 7^ = (rg_ £n 
(2) (2) 

I E q_, Ei ) /T E q_> Ei ■ Now inverting the external current 



I, T^ 3 ) changes sign and we get p' = P' E / P' E = (1 — 
73)/ (1 + 1q — 73)- From the above relations one can 

(3) (3) 

solve TV/ „ and T\,' „ provided that the probabil- 

(2) (2) 

ities, Tint and are known. One can evaluate T\ 
by applying Eq. JJJ) or the rates can be determined ex- 
perimentally. From Eqs. JSJ) and JJJ) we see that when 
a = (3, r^ 3 ) vanishes but remains finite. By keeping 

(2) 

AE fixed but setting B z = it is possible to measure T ■ 
independently. Figure |3| shows the asymmetry p — p' in 
the change of polarization with respect to the current in 
the source as a function of the magnitude of the current 
(bias voltage V = RI). 

In conclusion, we have studied the transitions induced 
by the third cumulant of current fluctuations on a probe 
quantum system. We have calculated a general formula 
for the transition rates and propose a scheme to measure 
the predicted results using a Josephson flux qubit. We 
have shown that the third-order transition rates are gov- 
erned by the variant of the third cumulant which to our 
knowledge has not been studied before. 
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